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Abstract 

This paper discusses two distinct, but related issues in quantum fluctuation effects. The 
first is the frequency spectrum which can be assigned to one loop quantum processes. The 
formal spectrum is a flat one, but the finite quantum effects can be associated with a rapidly 
oscillating spectrum, as in the case of the Casimir effect. The leads to the speculation 
that one might be able to dramatically change the final answer by upsetting the delicate 
cancellation which usually occurs. The second issue is the probability distribution for 
quantum fluctuations. It is well known that quantities which are linear in a free quantum 
field have a Gaussian distribution. Here it will be argued that quadratic quantities, such 
as the quantum stress tensor, must have a skewed distribution. Some possible implications 
of this result for inflationary cosmology will be discussed. In particular, this might be a 
source for non-Gaussianity. 
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I. INTRODUCTION 



This paper will deal with two aspects of quantum fluctuations. One will be the 
frequency spectrum which is associated with Casimir energy, that is, with the expec- 
tation value of the stress tensor operator or other quadratic operators. Some earlier 
work will be reviewed, in which it was shown that the frequency spectra are wildly 
oscillating functions, which nonetheless have finite integrals. The finite Casimir en- 
ergy always corresponds to a very small fraction of the area under one oscillation 
peak. The possibility of altering this remarkable cancellation will be discussed. 

The second aspect of quantum fluctuations to be considered is the probability 
distribution associated with the fluctuations of smeared operators. The well-known 
result of a Gaussian distribution for linear operators will be rederived. It will then 
be shown by way of a simple example that the fluctuations of quadratic operators 
can be described by a skewed, and hence non-Gaussian, distribution. Some further 
comments and speculations on the generic form of the probability distribution for 
stress tensor fluctuations will be offered. 



II. FREQUENCY SPECTRA 

It is well-known that the formal frequency spectrum of quantum fluctuations is 
flat; all modes appear on an equal basis in the expansion of a quantum field operator. 
This leads to the formally divergent zero point energy 

£o = E^- (!) 

However, this "unprocessed" spectrum is unobservable. The only quantities which can 
be observed arise from "processed" spectra, which have been modified by a physical 
process. The formation of a black hole, for example, processes the flat spectrum of 
incoming vacuum fluctuations and converts it into a Planck spectrum of outgoing 
particles, the Hawking radiation. The finite effects of one loop quantum processes 
may also be associated with a nontrivial frequency spectrum. Consider the Casimir 
effect as an example. The presence of boundaries modifies the divergent vacuum 
energy, Eq. by a finite amount, the Casimir energy. For a scalar field which 
satisfies periodic boundary conditions in one of three spatial dimensions, this energy 
is (Units in which h = c = 1 will be used in the remainder of this paper.) 

F - ** A (2) 

where L is the periodicity length and A is the transverse area. 

It is of interest to ask whether one can assign a finite frequency spectrum to the 
Casimir energy, that is, find a function ct(uj) such that 

E c = / duaiuj) . (3) 
Jo 
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This was done in Refs. where it was shown that one can obtain (t{lo) as a Fourier 

transform of the renormalized energy density operator at time-separated points. A 
more recent discussion is given in Ref. [3]. Let T^ v (t — t') be the renormalized stress 
tensor operator evaluated at two points separated in time by t — t'. Then 7)^(0), 
the operator evaluated at coincident points, is the observable stress-energy associated 
with the Casimir effect. In particular, 

E c = ALT a (0). (4) 

We can write 

^ ^ 1 Jut 



*( U ) = —I e^dt. (5) 



This function was evaluated explicitly for the electromagnetic case in Ref. |2|, and 
for the case of a scalar field in Ref. [lj. In the latter case, the result is 

a{ U ) = S{u,L) , (6) 

7T L 



where 

Six) = ]T (7) 



x3 sinfnx) 



n=l 71 



is a discontinuous, periodic function given by S(x + 2n) = S(x) and 

S( x ) = ^(tt-x), 0<x<2tt. (8) 

At first sight, the integral of cr(u) over all frequencies is poorly defined. However, 
it can be defined with a suitable convergence factor, such as an exponential function. 
One finds that the Casimir energy arises in the limit in which the convergence factor 
is removed. For example, 



oc 



E c = lim / du e'^ a(cu) . (9) 

/3^0 Jo 

The result is that the contributions of different frequency intervals almost exactly 
cancel one another, leaving a finite result which is small compared to the area of each 
of the peaks in Fig. [TJ 

A similar spectrum can be associated with the asymptotic Casimir- Polder potential 
between a polarizable particle, such as an atom in its ground state, and a perfectly 
reflecting wall. In the limit that the atom is far from the plate, compared to the 
wavelength associated with the transition between the ground state and first excited 
state, Casimir and Polder [^] showed that the interaction energy is, in Gaussian units, 

Vcp = ~P^, (10) 
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FIG. 1: The frequency spectrum, ct(uj), given by Eq. (jSJ) for the Casimir energy. The 
oscillations almost exactly cancel, leaving a net area under the curve equal to that of the 
shaded region indicated by the arrow. Here a is in units of 2A/L and uj in units of 2ir/L. 

where z is the distance to the wall, and cto is the static polarizability of the atom. It 
may be expressed as 



As before, the integral on uj may be denned using a convergence factor. The frequency 
spectrum in this case is illustrated in Fig. |21 

As we have seen, the frequency spectra associated with Casimir-type effects can 
be wildly oscillatory, and yet tend to integrate to a small net energy. An obvious 
question is whether it may be possible to slightly modify the contribution of specific 
frequency intervals and upset the delicate cancellation between the different positive 
and negative peaks in Figs. ^ and EJ For example, if one could alter the reflectivity of 
parallel plates just in a selected frequency range, it would seem to be possible to have 
a Casimir force much larger than that between perfectly reflecting plates, and which 
could be repulsive as well as attractive jij. A large, repulsive Casimir force could 
have significant industrial applications, such as allowing nearly frictionless bearings. 
Unfortunately, it is not at all straightforward to enhance Casimir forces. One might 
replace perfectly reflecting plates by dielectric slabs, which will have a finite, frequency 
dependent reflectivity. However, the Lifshitz theory {(J predicts that the force between 
two dielectric half-spaces will always be attractive and smaller in magnitude that the 
force in the perfectly reflecting limit. At least part of the reason for this is that the 





where, in this case, 



a(uj) = (2co 2 z 2 — 1) sin2cc;z + 2luzcos2luz . 
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FIG. 2: The frequency spectrum, cr(u), for the Casimir- Polder potential. The oscillations 
again almost exactly cancel, and the net area is equal to that of the shaded region indicated 
by the arrow. 

dielectric function s(cu) must satisfy the Kramers- Kronig relations, which follow from 
the requirement of analyticity in the upper-half u plane. This analyticity property is 
in turn a consequence of causality Q. Analyticity constrains e(u) in such a way that 
the contributions of various frequency regions continue to cancel efficiently. If one 
were to allow non-analytic dielectric functions, then it becomes possible to construct 
situations in which the Lifshitz thoery predicts repulsive forces 

Another approach to attempt to enhance Casimir forces was taken in Ref. 
There the force due to the quantized electromagnetic field on a small dielectric sphere 
near a perfectly reflecting plate was calculated A quasi-oscillatory result was obtained, 
which can be either attractive or repulsive, and is much larger in magnitude than the 
Casimir-Polder force, from Eq. pOjl . for a non-dispersive sphere. In Ref. it was 
shown that the qualitative form of the result does not depend upon whether the plate 
is perfectly reflecting or not. In both cases, the sphere was taken to have a plasma 
model dielectric function, and the resulting force oscillations as a function of distance 
from the plate have a scale set by the associated plasma wavelength. One can think of 
the effect of the sphere's nontrivial frequency response as upsetting the cancellations 
that occur for a non- dispersive sphere, and which are illustrated in Fig. |21 

If the vacuum modes of the electromagnetic field gave the only contribution to 
the force on the sphere, then there would be a large repulsive force on the sphere at 
certain separations which would be large enough to levitate the sphere in the earth's 
gravitational field, and should be observable. However, there is another contribution 
to the net force coming from quantum mechanical fluctuations of the electric charge 
within the sphere. This is the effect of what is sometimes called the plasmonic modes. 
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It was suggested by Barton [1(J| that the effect of the plasmonic modes should cancel 
the quasi-oscillatory terms coming from the vacuum modes. This suggestion is in fact 
correct, and was recently confirmed by an explicit calculation JjJ. There is still a 
nonzero net force, but it is always attractive and of the order of magnitude of that for 
a non-dispersive sphere. The cancellation occurs only if the sphere plasmon modes 
are in their ground state. If they are excited, then the quasi-oscillatory behavior 
appears in the net force. This is similar to the case of an atom. If the atom is in 
its ground state, then the force is attractive at all separations If the atom is 
in an excited state, however, the net force i s q uasi-oscillatory and typically much 
larger in magnitude than in the ground state |l2j]. It is still unclear why the plasmon 
modes can cancel the effects of the vacuum modes so efficiently, and whether one can 
engineer materials with repulsive Casimir forces without exciting plasma oscillations. 



III. PROBABILITY DISTRIBUTIONS 

In this section, we will turn to a somewhat different measure of vacuum fluctu- 
ations, the probability distribution associated with the fluctuations. First consider 
the case of a quantity which is linear in a free quantum field. It is well known that 
such quantities have a Gaussian distribution of probabilities. However, it will be 
worthwhile illustrating this result with an explicit example. Let <f(x) be a free scalar 
field and let 

<p = J if( x )dV (13) 

be a smeared field operator averaged over some spacetime region. The dimension of 
the spacetime is not relevant here. We could also include a sampling function f(x) 
and write f(x) dV in place of dV without changing the result. Next consider the 
expectation value of a power of (p in the vacuum state. All of the odd powers have 
vanishing expectation value because the n-point functions vanish for all odd n. That 
is, 

(f n+1 ) = J ((p 1<P2 ■ ■ ■ ^ 2n+1 ) dVxdV 2 ■ ■ ■ dV 2n+1 = , (14) 

where (pi = (p(xi), ect. However, all of the even moments are nonzero and can be 
expressed as powers of the second moment, 

(^ 2 ) = J . (15) 

For example, the identity 

{<Pi<fi<Pz<Pi) = {fif2) (^3^4) + (^1^3) (^2^4) + {fifi) (^2^3) (16) 
leads to the result 

<^ 4 > =3(^ 2 ) 2 . (17) 
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This is a special case of the general result 

=(2n-l)\\ (<?)». (18) 

The latter result can be obtained from the following counting argument. Wick's 
theorem allows us to decompose the 2n-point function (<pi<p2 ■ ■ • y?2n) into a sum of 
products of two-point functions, just as is illustrated in Eq. ([16)1 . All that we need 
to know is the number of terms in this sum. The first contraction of the 2n-point 
function contain 2n — 1 terms, as we can select any field to start, and it then has 
2n — 1 partners with which it can be contracted. Similarly, the contraction of each 
(2n — 2)-point function to (2n — 4)-point functions contains 2n — 3 terms, and so on, 
leading to the factor of (2n — 1)!! in Eq. (|18j) . 

This result for the general even moment implies that the probability distribution 
is a Gaussian function 

P{<p) = — L= e-^/^> , (19) 
which is confirmed by the facts that 

/oo 
x 2 P(x) dx (20) 
-oo 

and that 

poo 

{(p 2n )= x 2n P(x)dx = (2n - l)\\ (<p 2 ) n . (21) 



The result that P((p) is a Gaussian function is responsible for the prediction of Gaus- 
sian fluctuations in inflationary cosmology. In most versions of inflation, density 
perturbations are linked to the quantum fluctuations of a scalar inflaton field, which 
is treated as a free field [Tj3]. This leads to a Gaussian distribution of density fluctu- 
ations, which seems to be consistent with observation |l4j . 

Let us now turn to the fluctuations of quadratic operators, such as the smeared 
stress tensor. Now there is no particular reason to expect the probability distribution 
to be symmetric, much less Gaussian. In fact, as we will see from an explicit example, 
the distribution is in general a skewed one. The example will involve a massive scalar 
field in two-dimensional Minkowski spacetime, but in the limit of a very small mass. 
The Hadamard function for this field is 

1 1 



G(x, x') = - {<p(x), y{x')} = -- N (mVAt* - Ax 2 ) , (22) 

where m is the mass and N Q is a Neumann function. In the limit of small m, this 
function becomes 

G(x, x') = ~ ln[/i 2 (At 2 - Ax 2 )] , (23) 

4:71 

where /i = e 7 m/2 and 7 is Euler's constant. The commutator function, in the limit 
that m — > 0, is 



G c {x, x') = Mx), cp{x')) = % - [6{Ax - At) - 9(Ax - At)} , (24) 
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where 9{x) = 1 for x > and 8(x) = — 1 for x < 0. Note that the commutator 
function is finite for m — 0, whereas the Hadamard function has a logarithmic di- 
vergence in this limit. This is a well-known feature of the massless scalar field in 
two-dimensional Minkowski spacetime. This divergence can be removed by selecting 
a vacuum state which breaks Lorentz invariance |l5[ , but here we will assume a small, 
nonzero mass. 

Next consider the average of the normal ordered square of (p over a finite time 
interval 



T Jo 

In this model, all of the field operators will be taken to be at the same point in space, 
so that Ax = 0. We wish to calculate the vacuum expectation value of various powers 
of (p 2 . Begin with the second moment, 

<(^ 2 ) 2 > = Y2 [ ^ dt ^ ^? : : ^ : > = h la dh lo dt2{lf{tl) V{t2))2 ' 

(26) 

where the last step follows from Wick's theorem. The unsymmetrized two-point 
function can be written as a sum of the Hadamard and commutator functions 

{tpfa) <f(h)) = G(ti, t 2 ) + X - G c (h,t 2 ) . (27) 

We can now evaluate the second moment explicitly and find 

(M 2 > = ^ ln 2 (T/,) (28) 

in the limit that Tfi <C 1. In this limit, the Hadamard function part of the two-point 
function gives the leading contribution. 
Next we turn to the third moment 

<(^)3) = _L f T dtl F dt2 F dt 3 (: ifit.f : : V {t 2 f :: <p(t 3 ) 2 :) , (29) 
I 6 Jo Jo Jo 

which can be written as 

((^) 3 ) = — / dtj dt 2 f ^ 3 M*l)^2))W 2 )^3))Mil)^3))- (30) 

1 Jo Jo Jo 
In the limit of small T/i, this becomes 

<(^ 2 ) 3 > = ~ ln3 ™ 

7T° 



7T VT« 



(31) 



The third moment is positive ((yj 2 ) 3 ) > 0, and of the same order of magnitude as the 
second moment in the sense that 

((^) 3 )^ = v / 2((^ 2 ) 2 )^. (32) 
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Thus the probability distribution is significantly skewed. 

The corresponding calculations for components of the quantum stress tensor have 
not yet been performed, but it is reasonable to guess that the results will be similar 
to those of this ip 2 model in two dimensions. The technical details are somewhat 
more complicated, however. Let 

f = J:T tt : f(x)dV (33) 

be the average of the normal-ordered energy density over a spacetime region defined 
by the non-negative sampling function f(x). The stress tensor correlation function 
is singular in the limit of null separated points, 

(:T u (x)::T u (x'):)~(x-x')- 2d , (34) 

where d is the number of spacetime dimensions. Thus the integral for the second 
moment (T 2 ) will appear to contain a non-integrable singularity. However, it is 
possible to define such integrals using an integration by parts procedure 3, 12 > 3- 
It is also possible to use dimensional regularization [19(. In both approaches, the 
integral becomes finite. The third and higher moments can be defined in a similar 
manner. If (T 3 ) > 0, as in the case of the ip 2 mode, then the probability distribution 
will again be skewed. There should be a connection between this distribution and the 
quantum inequalities which set lower bounds on the expectation value of T in non- 
vacuum states ji^, H3, 0, 23L 24l 25l 26l 27. 128^ . In the case of temporal averaging 



only, these inequalities are of the form 



(T>>Tn, = — , (35) 



where r is the characteristic width of the sampling function and C is a constant. 
Thus negative energy densities are allowed in quantum field theory, but are tightly 
constrained. The lower bound in the expectation value in an arbitrary quantum 
state, and the lower bound on the probability distribution of fluctuations in the 
vacuum state are expected to coincide |29j |. One way to understand this is to imagine 
writing the Minkowski vacuum state as a superposition of an alternative basis set of 
states. In each of these states, the measured T is bounded below by Tib- Hence the 
possible fluctuations which can be observed in the Minkowski vacuum state are also 
bounded below by the same value. If the probability distribution P(T) has positive 
skewness and a finite negative lower bound at T = Tib, then we can make some general 
observations on P(T), The first moment is defined to vanish, (T) = 0. A positive 
third moment implies a long tail in the positive direction. Thus the greater portion 
of the area will lie to the left of T = 0. Thus a typical measurement of the sampled 
energy density is more likely to yield a negative than a positive value. However, when 
the measured value is positive, its magnitude is likely to be greater than when it is 
negative. Figure|3]is a sketch of a hypothetical probability distribution which satisfies 
all of these conditions. A more detailed analysis of the probability distribution for 
various cases is currently in progress ji^j. 
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FIG. 3: A possible probability distribution P(T) for stress tensor fluctuations is illustrated 
as a function of an averaged stress tensor component T. This distribution is zero below the 
lower bound, T < Tib < 0, but has a nonzero tail extending infinitely far in the positive 
direction. The majority of the area under the curve lies in the negative region, T < 0. 

The non-Gaussian nature of quantum stress tensor fluctuations may have applica- 
tions to inflationary cosmology in the form of a source of a non-Gaussian component 
in the density perturbations. This is currently under investigation |30(. 
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